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O . Abstract 

The gluon propagator is investigated in Landau and in maximal 
center gauge for the gauge group SU(2) by means of lattice gauge sim- 
^j^. ulations. We find that Gribov ambiguities arising from the implemen- 

Q^l tation of Landau gauge have a small influence on the propagator. In 

^ ' agreement with previous findings, we obtain that the small momentum 

behavior is dominated by a mass, which is of order (1.48 it 0.05) -^/a, 
where a is the string tension. By removing the confining vortices from 
^ ' the full Yang-Mills ensemble, we convert full YM-theory into a the- 

j^ ■ ory which does not confine quarks. We find that in the latter case 

the strength of the gluon propagator in the intermediate momentum 
range is strongly reduced. The spectral functions which reproduce the 
numerical data for the propagators are analyzed using a generalized 
Maximum Entropy Method. 
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1 Introduction 

Nowadays, the general belief that QCD is the correct theory of strong in- 
teractions is supported by numerous high energy collisions experiments (see 
e.g. [0]). At high momentum transfers, at which the effective quark gluon 
coupling becomes small, predictions of perturbation theory nicely agree with 
the experimental data. On the other hand, this coupling becomes large and 
even diverges at low energies (Landau pole) thus prohibiting a perturbative 
treatment and pointing towards new phenomena. One might speculate that 
the strong increase of the effective coupling is the progenitor of quark con- 
finement. Note, however, that a toy quark model which mimics the QCD 
perturbative behavior of the running coupling strength reveals that the new 
phenomenon at low energies is quark mass generation (which screens the 
Landau pole P|) rather than quark confinement. 

One prominent method to treat non-perturbative Yang-Mills theory is the 
numerical simulation of lattice gauge theory (LGT). LGT covers all non- 
perturbative effects and, in particular, bears witness of quark confinement 
(see e.g. [||). Moreover, the numerical results suggest that, in certain gauges, 
topological defects act as quark confiners. The presence of these defects 
is independent of the coupling strength, and their properties are beyond 
the reach of perturbation theory. In the case of Abelian gauges [||, these 
defects are chromomagnetic monopoles, which condense and give rise to the 
dual Meissner effect as confinement scenario p, H, ^ (for a caveat see 0). 
A generalization of the dual Meissner effect to the non-abelian case was 
put forward in |^ to resolve the neutral particle problem. In the Center 
gauges @], the topological defects are center vortices. The mechanism of 
quark confinement can be understood as a percolation of center vortices 
which acquire physical relevance in the continuum limit 0. Lattice gauge 
calculations also provide an intuitive picture in terms of vortex physics for the 
deconfinement phase transition at finite temperatures |jlO|]. Reducing the full 



Yang-Mills configurations to their vortex content still yields the full string 
tension [§]. Vice versa, removing these vortices from the Yang-Mills ensemble 
results in a vanishing string tension P] and in addition to a restoration of 



chiral symmetry |11 



Unfortunately, simulations of LGT including dynamical quarks are still cum- 



bersome despite the recent successes by improved algorithms |]T2[ and the 
increase of computational power. Moreover, at the present stage systems at 
finite baryon densities are hardly accessible in the realistic case of an SU(3) 
gauge group |]T3[ (for recent successes see W^). These deficiencies are over- 



come by a second non-perturbative method to treat Yang-Mills theory: the 
approach by the Dyson- Schwinger equations (DSE). By contrast to the lattice 
formulation, the DSE approach can easily deal with dynamical quarks and, 
furthermore, can be easily extended to finite baryon densities |TH] . Moreover, 



it can be used to study hadron phenomenology |I6|, 0. The disadvantage 
of the DSE approach is that it requires a truncation of the infinite tower of 
equations, and that this approximation is difficult to control and to improve 
systematically. In addition, the DSE approach needs gauge fixing which is 
obscured by Gribov copies. Whether the standard Faddeev-Popov method 
of gauge fixing is appropriate in non-perturbative studies, is still under de- 
bate 



In view of the importance of the DSE approach for the understanding of 
hadron physics []T^ |T^, it is highly interesting to learn more about the im- 



pact of truncating the Dyson tower of equations. Since the gluon and ghost 
propagators are essential ingredients of the quark DSE, these propagators 
are of particular interest. 

In this paper, we will address the gluon propagator in Landau gauge for the 
simpler case of pure SU(2) lattice gauge theory. The lattice result will be 
compared with the one provided by the solution of the coupled ghost-gluon 



Dyson equation |T9|, ^, |21], ^. By removing the confining center vortices 
from the SU(2) Yang-Mills ensemble, we will focus on the information on 
quark confinement which might be encoded in the gluon propagator. High 
accuracy numerical data for the latter are obtained by means of a new nu- 
merical method superior to existing techniques. 

The organization of the paper is as follows: In section 2, we present the 
lattice definition of the gluon propagator. For making contact with the ab 
initio continuum formulation of Yang-Mills theory, we will define the gluon 
field from the adjoint links. We furthermore outline in this section the gauge 
fixing procedure. In section 3, we present our numerical results for the gluon 
propagator in the Landau gauge and the maximal center gauge, respectively. 
In section 4, we discuss the gluonic spectral function which is extracted from 
the gluon propagator by using (an extension of) the maximal entropy method. 
Finally, our conclusions are presented in section 5. 



2 The lattice approach to the gluon propa- 
gator 

In this section, we will extract the gluon propagator of the continuum Yang- 
Mills theory by considering the continuum limit of the lattice gauge theory. 
Thereby, we will carefully examine the relation between the lattice link vari- 
ables and the gauge potential. 

2.1 The lattice definition of the gluon field 

Before identifying the gluonic degrees of freedom in the lattice formulation, 
we briefly recall the definition of the gluon field in continuum Yang-Mills 
theory. For simplicity, we will consider the case of SU(2) Yang- Mills theory. 

Under a gauge transformation of the fundamental matter field, i.e. 

q{x) -^ q'{x) = n{x) q{x) , n{x) G Sf/(2) , (1) 

the gluon fields transforms as 

A^^'ix) = O^^x) A'l^ix) + -t''^^ 0'"'{x)d^Of\ (2) 

0''\x) := 2tr{fi(s)t"fit(x)t^} , 0''\x) E ^0(?,) . (3) 

Let us stress that the gluon fields transform according to the adjoint represen- 
tation while the matter fields are defined in the fundamental representation. 

Let us compare these definitions of fields with the ones in LGT. In LGT, 
a discretization of space-time with a lattice spacing a is instrumental. The 
'actors' of the theory are SU(2) matrices U^{x) which are associated with the 
links of the lattice. These link variables transform under gauge transforma- 
tions as 

Ul{x) = Vl{x)U^{x)n\x + ^i) fi(a;) G Sf/(2) . (4) 

For comparison with the ab initio continuum formulation, we also introduce 
the adjoint links 

U;\x) := 2iv{U,{x)eul{x)t'], t'^ = V (5) 

U'^{x) := 0(x)W^(x)0^(x + /i) , 0(x)gS0(3), (6) 

where 0{x) was defined in (H), and r" are the Pauli matrices. 



In order to define the gluonic fields from lattice configurations, we exploit 
the behavior of the (continuum) gluon fields under gauge transformations 
(see (0)), and identify the lattice gluon fields v4^(x) as algebra valued fields 
of the adjoint representation, i.e. 

W^V) =: [exp{t^4(x)a}J , tl := e'^^% (7) 

where the total anti-symmetric tensor e"*'^ is the generator of the SU(2) group 
in the adjoint representation, and where a denotes the lattice spacing. 

For later use, it is convenient to have an explicit formula for the (lattice) 
gluon fields A?^(x) defined by (0) in terms of the SU(2) link variables U^{x). 
Usually, these links are given in terms of four-vectors of unit length 

U,{x) = u^ix) + zu,{x)t, [u'^ixf + [u,{x)f = 1 . (8) 

Inserting this representation for U^{x) into (|), we expand (|) and (0), re- 
spectively, in powers of the lattice spacing using 

M^(x) = 0(a), iuiyix) = l-Oia'), 

Comparing the order 0{a), we find 

A^^{x) a + 0(0^) = 2M°(a;) m|^(x) , without summation over /i . (9) 

Let us emphasize that this gauge field v4|^(x) has been defined by the adjoint 
link (^. As a consequence, the representation (P) is invariant under a non- 
trivial Z2 center transformation, i.e. U^{x) -^ — f/^(x). 

Let us contrast our definition (|) with the previous definition of the gauge 
field from the fundamental link, i.e. 

U^{x) = exphaA^t^l , aA^ = 2uJ.(x) + ^(a^), (10) 

where one must assume that the link field U^{x) is close to the unit element, 
i.e. u°(x) = 1 — Oi^o?). Indeed, the gluon field A^^{x) changes sign under 
a non-trivial center transformation. The previous definition of the gauge 
field A^,, therefore contains information on center elements and coset fields 
as well. Here, we propose to disentangle the information carried by center 
elements and coset fields A^^{x) and to study their correlations separately. In 
the following, we will present the correlation function of the coset "gluon" 
fields Al{x). 



2.2 Gauge fixing 

Calculation of the gluon propagator requires gauge fixing. In order to be able 
to compare with the Dyson- Schwinger approach, we firstly use the (lattice) 
Landau gauge condition 

n{x) : ^ trf/;,(x) -> max, (11) 

{x},fj, 

where U' [x) is the gauged transformed link (^. Gauge transformation Q{x) 
which are determined by (0) bring the SU(2) link elements as close as pos- 
sible to the unit element. Decomposing the link variable in this gauge as 

U'^ix) = Z,{x) exp{^4(x)t^a} , (12) 

where 



Z,{x) = signtrf/^(x) G {-1,+1} , JE^m^m ^ [0,-] , (13) 



one observes that in (lattice) Landau gauge (|TT|) the role of the Z2 center 
elements are de-emphasized (almost all Z^j_{x) are 1) and most of the physics 
is contained in the adjoint field A^^{x) (^). For this reason, we do not expect 
a vastly different gluon propagator when the more standard definition of the 
lattice gluon fields, i.e. A^^{x) (see (0)), in terms of the fundamental links 



U^{x) is used [|3], H 



In the case of the gauge (|Tl|) , the gluon fields (H) satisfy the familiar Landau 
gaugeQ 

d,A^\x) = . , (14) 

It is this gauge which is used in the Dyson-Schwinger studies. In this ap- 
proach, the gauge condition ( [T^ ) is implemented by means of the Faddeev- 
Popov method, thereby relying on the assumption that the Faddeev-Popov 
determinant represents the probabilistic weight of the gauge obit specified 
by its representative A'°'{x). This method is correct if the gauge condition 
picks a unique solution VL{x) of ([l^) for a given field A^^{x). Unfortunately, 
the Landau gauge condition generically admits several solutions depending 
on the "background field" A\^{x) (Gribov ambiguity). Further restrictions on 
the space of possible solutions VL{x) are required |^. It was argued in |TS 



Strictly speaking, the eq.([14[) holds up to singular points where Z^{x) = — 1 (see 
. These singularities are subject of future work |23]. 



that the Faddeev- Popov method is not always justified if Gribov ambiguities 
are present. 

Let us contrast the continuum gauge fixing with its lattice analog. In a first 
step, link configurations Uf^{x) are generated by means of the gauge invariant 
action without any bias to a gauge condition. In a second step, the gauge- 
fixed ensemble is obtained by adjusting the gauge matrices Q{x) (see (|)) 
until the gauged link ensemble satisfies the gauge condition (p!T|). When 
one representative of each gauge orbit is picked by this procedure, one must 
determine with which weight this representative contributes to the observable 
of interest. In the absence of Gribov ambiguities, this weight factor is given 
by the Faddeev Popov determinant (see e.g. |^). The important point is 
that in the lattice calculation the Faddeev Popov determinant needs not to be 
explicitly evaluated. This is because the unbiased generation of gauge field 
configurations which are subsequently transformed into the desired gauge 
produces each configuration of the gauge orbit with equal weight. Hence, 
this procedure automatically produces the proper probability distribution of 
the gauge fixed sub-manifold to which the representative belongs. Further 
details of the numerical approach to gauge fixing can be found in appendix 

E 

However, this does not dispense us from dealing with the Gribov ambiguity. 
Let us illustrate this point: The naive Landau gauge condition for the gluon 
field (|l^) is satisfied if we seek an extremum (instead of the maximum) of 
the variational condition (0). If we restrict the variety of solutions VL{x) 
which extremize ( pH]) to those solutions which maximize the functional (pJ]), 
we confine ourselves to the case where the Faddeev-Popov matrix is posi- 
tive semi-definite. The corresponding fraction of the configuration space of 
gauge fixed fields A'^ is said to lie within the first Gribov horizon. How- 
ever, there is still a variety of possible gauge transformations VL[x) which all 
correspond to local maxima of the functional (|l3)- A conceptual solution 
which resolves this residual Gribov ambiguity is to restrict the configuration 
space of gauge fixed fields A^ to the so-called fundamental modular region. 
In the lattice simulation, this amounts to picking the global maximum of 
the variational condition (pTI). In practice, finding the global maximum is 
a highly non-trivial task. A numerical algorithm which obtains the gauge 
transformation matrices VL{x) from the condition (^) might fail to locate 
the global maximum, and the numerical simulation might still sample a par- 
ticular set of local maxima. Different algorithms might lead to different local 
maxima, and, hence, implement different gauges. A comprehensive study of 
the Gribov problem using lattice methods can be found in [^ . 



Here, we will study two extreme cases of gauge fixing: firstly, we will im- 
plement a gauge by means of an iteration over-relaxation algorithm which 
almost randomly averages over the local maxima of the variational condi- 
tion ([TT|) (10 gauge). This is the standard gauge fixing algorithm commonly 



used by the community (for technical details, see e.g. ||2^). The resulting 
gluon propagator will then be compared with the gluon propagator of a gauge 
where a simulated annealing algorithm searches for the global maximum, (SA 



gauge). In the latter case, we used the algorithm outlined in [^. We stress 
that it is not granted that the SA algorithm finds the global maximum, i.e. 
that the gauge is fixed to the fundamental modular region. Instead, another 
set of local maxima is randomly chosen. The SA approach generically yields 
higher values of the gauge fixing functional than the 10 algorithm. It is the 
present state of the art to study the degree with which the observable of 
interest depends on the choice of the these two extreme cases of gauge fixing. 
This will be done below. It will turn out that the gluon propagators of both 
gauges agree within statistical error bars. These observations indicate that 
the gluon propagator does not depend on the subset of configurations which 
we choose from the first Gribov domain. 

2.3 Form factor calculations 

The link configurations are generated using the Wilson action. We refrain 
from using a "perfect action" since we are interested in the gluon propagator 
in the full momentum range; simulations using perfect actions recover a good 
deal of continuum physics at finite values of the lattice spacing at the cost of 
a non-local action. For practical simulations, perfect actions are truncated 
which becomes an un-justified approximation at high energies where the full 
non-locality of the action must come into play. 

In the present paper, most of the calculations were performed using a 16^ x 32 
lattice. The dependence of the lattice spacing on /? (renormalization) , i.e. 

aa2(/3) = 0.12 exp|-—(/3- 2.3)}, a := {AAOMeVY , (15) 

is appropriate for j3 G [2.1,2.6] for the achieved numerical accuracy. Given 
([T5|) , it is straightforward to calculate the extension of the lattice in one 
direction, i.e L^ = N^a{l3), where N^ is the number of lattice points in x- 
direction. In order to estimate the momentum range which is covered by the 
actual simulation, it is convenient to introduce the UV-cutoff by A = 7r/a(/5). 
Table |l| provides Lt and A for Nt = 32. 



Once gauge-fixed ensembles are obtained by implementing a variational gauge 
condition (see discussion of previous section), the gluon propagator is calcu- 
lated using 

D;lix-y) = {A';^{x)Al{y)),,c, (16) 

where ^J^(a;) is defined in terms of the coset part of the link (see (^). The 
Monte-Carlo average is taken over 200 properly thermalized gauge configu- 
rations. Of particular interest is the Fourier transform of this propagator 
which is defined by 

DtM = «' E ^tui^) expjzpa:} , p, := j^^n, , (17) 

X 

where Uk labels the Matsubara mode in fc-direction and where A^^ is the 
number of lattice points in this direction. It is also convenient to define the 
lattice momentum p by 

which coincides with the Matsubara momentum p^ in (0) in the limit n^ <^ 
Nk- The definition (|TBp has the advantage that e.g. the free lattice propagator 
takes the familiar form 1/pkPk- 

From perturbative Yang-Mills theory, one expects that the deviation of the 
full propagator from the free one is logarithmically small for large p"^. In 
order to work out the non-trivial information of the Yang-Mills dynamics on 
the propagator, it is useful to introduce the gluon form factor F(p^) by 



D 

p 



^, D{p) := J2d;1{p) (19) 



a,fj, 



which measures the deviation of the full propagator from the free one. Since 
in Landau gauge the propagator is diagonal in color and transversal in 
Lorentz space, the form factor F{p^) contains the full information. 

2.4 Numerical method 

In principle, the gluonic form factor can be obtained by the Fourier trans- 
form of the measured quantity (|16D. Note, however, that the propagator ([T6|) 
is a rapidly decreasing function in coordinate space implying that the phys- 
ical information at large distance |x — ?/| is washed out by statistical noise 
and that, therefore, the information on the low momentum behavior is lost. 

9 



Table 1: Simulation parameters 



/? 


2.1 


2.2 


2.3 


2.4 


2.5 


L [fm] 


8.6 


6.6 


5.0 


3.8 


2.9 


A [GeV 


2.3 


3.0 


4.0 


5.2 


6.8 



To overcome this deficiency, it was proposed in |3^ to directly address the 
propagator in momentum space, i.e. 



1 



^(^) = J2j;p\ [5Z^m(^) cos(px)J + [J2^Uy) ""'^(py) 

r, II Siii ^ ^1 



a,fj, 



MC 



(20) 

where A'^aii is the number of lattice points. By the help of the translational 
invariance, i.e. 

{A'^^{x)Aliy))^^ oc f{x-y), 

one indeed finds from (Wi) that 



D{p) = E E ( ^m(^) ^m(o) ;^,^ ^o^ I p^ 



(21) 



a, II X 



Since the function D{p) also contains the trivial factor 1/p^, which is also 
present in a free theory, a further increase of the numerical accuracy is 
achieved by directly addressing the form factor F{p'^). For these purposes, 
we firstly choose the momentum transfer to p = (0, 0, 0,p4)^ without any loss 
of generality and define 



AiA^(x) := A^{x + aei) - A^{x) 



:221 



where 64 is the unit vector in time direction. A straightforward calculation 
yields 



^(^4) = E]^( [E^*^m(^)cos(^^) 

n 1 1 3iii ry 



a,fj, 



+ [j2^tA';{y)sm{py) 
y 



MC 

(23) 

where p^ is the lattice momentum (^) in time direction. By sake of the Aj 
operators in ( ^3] ) the free part, i.e. 1/p^, of the propagator D{p) is precisely 
canceled, and we are left with the quantity of interest F{p^). It turns out 
that Monte- Carlo average (p3D allows for a high precision measurement of 
the form factor. 
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Figure 1: The gluonic form factor -F(p^) as function of the momentum trans- 
fer (left panel: linear scale; right panel: log-log scale). Also shown is the 
solution of the set of DSEs proposed in 
case of SU{2) 



19 which have been solved for the 



3 Numerical Results 



3.1 Landau gauge 

In this section, we will mainly use the iteration over-relaxation algorithm 
to implement the Landau gauge condition (]Tl|), and will refer to the corre- 
sponding gauge as lO-gauge. Later, we will also compare this gluonic form 
factor with one obtained in the SA-gauge where the gauge condition ( ]TT| ) is 
implemented by a simulated annealing algorithm. 

In order to obtain the form factor -F(p^) as function of the (lattice) momen- 
tum p (p!8|), 200 properly thermalized configurations of a 16^ x 32 lattice 
were used. Physical units for the momentum p (|T8D can be obtained by 
using ([T5|) . Calculations with different /3- values correspond to simulations 
with a different UV-cutoff A := 7r/a(/3). In the first place, we obtain the un- 
renormalized form factor Fb(p^) as function of the momentum in physical 
units. The desired renormalized form factor Fr{p^) is obtained via mul- 
tiplicative renormalization, i.e. Fji(p'^) = Z^^Fb{p^)- Thereby, the gluonic 
wave function renormalization .^3 (A) is chosen to yield a finite (given) value 
for the renormalized form factor Fr at a fixed momentum transfer (renor- 
malization point). In the following, we will suppress the subscripts of the 



11 



form factors, and F{p'^) always refers to the renormalized form factor. 

Figure |I| shows the result for the gluonic form factor in the lO-Landau gauge. 
Our reference scale is the string tension, which we set to a = (440 MeV)^ in 
order to assign physical units to momenta. The data for the renormalized 
form factor -F(p^) which were obtained with different /3-values nicely agree 
within numerical accuracy, thus signaling independence from the UV-cutoff 
and the lattice volumes in the range shown in table ^ 

At high momentum the lattice data are consistent with the behavior known 
from perturbation theory. 



Fip") oc 1/ 



log 



p" 



2 1 13/22 



P^ ^ /i^ > (1 GeV)^ 



(24) 



Also shown in figure |l| is the coarse grained "mass fit" (A/", mi, ■m2, "hil, s 
fitted parameters; momentum p and all mass scales in units of 1 GeV) 



Fip" 



N- 



p^ 



p^ + mf 



p'* + 7712 



+ 



1 13/22 



(25) 



[log {ml +p'^) J 

which nicely reproduces the lattice data within the statistical error bars for 
A/" = 8.1133, mi = 0.64 GeV, mg = 1.31 GeV , (26) 



8.1133 
0.32, 



mi = 0.64 GeV 
rriL = 1.23 GeV . 



Note that the multiplicative renormalization only affects the normalization 
A/" implying that the other fitting parameters are renormalization group in- 
variant quantities. 



on 



Let us point out that an effective gluonic mass was introduced in p2 
phenomenological grounds. Gluonic masses were realized as electric and 
magnetic screening masses in the high temperature phase of SU(2) YM- 

and 



theory |3^, and were also reported in the maximal Abelian gauge 

in the Laplacian Landau gauge [^. Masses have been also reported for the 



case of a SU(3) gauge group in Landau gauge |gj and in Laplacian Landau 
gauge ||35|, ^ . Note, however, that the description of the propagator in terms 
of two conjugate mass poles seems adequate for the Coulomb gauge |]37 |. 



3.2 Comparison with DSE solutions 

Despite more than twenty years of successful hadron phenomenology origi- 
nating from the quark |T^, |T^ and gluon [^ ^ ^ DSE, respectively, the 
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coupled set of continuum DSEs for the renormalized gluon and ghost propaga- 



tors has only recently been addressed in [jT9[ and subsequently in [^, ^ |22| 



In ref. [|TS|, it was firstly pointed out that, at least for a specific truncation 
scheme, the gluon and ghost from factors satisfy scaling laws in the infra-red 
momentum range, in particular 

F{p')^[p'Y\ p'<^AlcD- (27) 

Depending on the truncation of the tower of Dyson equations and on the an- 



gular approximation of the momentum loop integral, one finds k = 0.92 |19] 



or K = 0.77 [^ or fv — >■ 1 |]2T|. The lattice data are consistent with k = 0.5 
corresponding to an infra-red screening by a gluonic mass (see figure |l| right 
panel). Interestingly, the prediction that the running coupling strength 

«(/) = a{ij'^) F(p2) G^(/) , G'(p^) : Ghost form factor (28) 

approaches a constant in the limit j9^ ^ is independent of the truncation 



and approximations used in [19, EOl UM, O 



A feature of full Yang-Mills theory is multiplicative renormalizability (MR) 
which guarantees that the form factors F{p^) and G{p^) might be rescaled 
independently to satisfy the renormalization condition F{ij^) = 1, Gi^fi"^) = 
1 (yU renormalization point). Above, we have made use of this feature to 
obtain the renormalized form factor FijP') from that lattice data for several 



values of (3. It turned out that the truncations studied so far in |19 



2l| violate multiplicative renormalizability to a certain extent (see p^ and 



22| for detailed discussions). Progress was made in [^ where a system of 
renormalized coupled gluon ghost DSEs were derived which manifestly reflect 
MR. Unfortunately, a self-consistent solution to this set of DSEs has not yet 
been obtained. 

For a quantitative comparison of our lattice data with the SU(2) DSE so- 
lution, we refer to the truncation scheme of [|Ty] which incorporates Taylor- 
Slavnov identities to some extent. In this approach, the running coupling 
at the renormalization point /x serves as an input to assign physical units to 
the momenta. Assuming an approximate MR, the solution for -F(p^) is then 
rescaled to satisfy -F(/i^) = 1. The final result is compared with the form 
factor obtained by our lattice calculations. 

The DSE approach requires a knowledge of the scale Aqcd which is (to one 
loop accuracy) defined by (see e.g. |42[| ) 



In(pVA^CD) 
13 



which holds for p^ ^ /x^ ^ ^%cd- Defining 

/i^ dF{p') 



P 



and using (pQ]), a straightforward cafoulation gives 

13 



(30) 



AlcD = -"' ^^P ^ ^ 1^ • (31) 



For the parameter set (|2b| ) which fits the lattice gluonic form factor with the 
string tension a = (440 MeV)^ as reference scale, p can be estimated from 
the fit function (l25|) . The renormalization point, fi = 6 GeV, was chosen to 
be part of the asymptotic momentum regime where both non-perturbative 
methods, i.e. LGT and the DSE approach, reproduce the known perturbative 
behavior. We finally obtain Aqcd ~ 889 ±10 MeV, which is of the same order 
of magnitude than its SU(3) analog when the SU(3) gauge theory is equipped 
with the same string tension. At the renormalization point /i = 6 GeV, the 
running coupling constant for the SU(2) gauge theory is obtained by 

Alt 

a(/i = 6GeV) = — , , , ,,.^ ^ 0.449 (32) 

^ ' 22/3 ln(/i7A|^^) ^ ^ 



Using (|3^ as input, the solution of the DSEs with the truncations of [19 



was solved [^ for the case of the SU(2) gauge group which yields the DSE 
result shown in figure |I|. We find a qualitative agreement of the DSE solution 
with our lattice result. Given the varieties of exponents k ( P?] ) depending 
on the truncations of the DSEs, we do not expect a detailed agreement of 
the form factor close to zero momentum. We point out that the peak of 
the form factor at the intermediate momentum range is also observed in 
the DSE approach, although there are quantitative deviations concerning its 
steepness. 

3.3 Gluon propagator and confinement 

In the so-called maximum center gauges |g|, the role of the coset part of 
the links referred to as "gluons" for the infra-red physics is de-emphasized, 
and center vortices appear as physical degrees of freedom in the continuum 
limit 0. In this gauge these vortices act as the confiners of the theory: 

^ We thank C. Fischer for communicating his DSE solution for the SU(2) case prior to 
pubhcation. 
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Figure 2: The static quark anti-quark potential (left panel) and the corre- 



sponding gluonic form factors (right panel); DSE solution from |]3T |. 



reducing the full Yang-Mills configurations to their vortex content still yields 
the full string tension [R|, and, vice versa, removing these vortices from the 



Yang-Mills ensemble results in a theory with a vanishing string tension 0, 11 



In order to get a handle on the information of quark confinement encoded 
in the gluon propagator in Landau gauge, we firstly modify the SU(2) Yang- 
Mills theory to a theory which does not confine quarks, and, secondly, we 
compare the form factor of the gluon propagator of the modified theory with 
the full SU(2) result (see figure |I|). 

Let us briefly outline the numerical procedure for removing the confining cen- 
ter vortices from the Yang-Mills ensemble: firstly, we implement the maximal 
center gauge (MCG) condition 



^x): $^[trf/;(. 

{x},/i 

by an iteration over-relaxation algorithm 
by center projection 



max 



(33) 



and obtain the 2' (2) vortex links 



V,{x) 



MCG JJMCG^^^ Pioj ^^^^^ 



sign trf/ 



MCG I 



The elements Z^{x) span an Z{2) gauge theory which contains the confining 
vortices as physical degrees of freedom. In order to remove the vortices from 
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Figure 3: The gluonic form factor -F(p^) in MCG as function of the momen- 
tum transfer (left panel: linear scale; right panel: log-log scale) for the case 
of lO-MCG. 



the full ensemble, we define configurations 



ur'^(x) 



Z^{x) 



U,{x) 



The configurations f/™ (x) do not yield quark confinement any more: the 
static quark anti-quark potential calculated from ensembles U'^°'^{x) (see 
figure 1^ left panel) illustrates that a removal of the center vortices produces 
a non-confining theory. 

Finally, we implement the lO-Landau gauge condition on the configurations 
UJ^°'^{x) and calculate the corresponding renormalized form factor by the 
procedure outlined in the previous section. Figure ^ right panel shows this 
gluonic form factor obtained from the modified ensemble. The striking fea- 
ture is that the strength of the form factor in the intermediate momentum 
range is drastically reduced. 



3.4 Maximal center gauge 



In the maximal center gauge |P|- [Tl|], important parts of the infra-red physics 
is shifted to the vortex degrees of freedom. It is therefore instructive to 
investigate the residual information carried by the gluonic form factor in this 
gauge (P^l). We used an standard iteration over- relaxation algorithm M to 



determine the gauge matrices f2(x) for a given "background" configuration 



U,{x) 
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Once the gauge condition (^) is implemented, we used the adjoint hnk to 
define the gluon field (see @). The resulting gluon field satisfies the Landau 
gauge condition d^^A^ (x) = 0, since the MCG is equivalent to the adjoint 
Landau gauge 

\r^ r 1 ^ 

Q{x) : 22 h^^M(^) ^ ^^^ ' (34) 

where W^{x) is the adjoint link (|]). 

A thorough study of the MCG in the continuum limit was performed in [^ . 
One finds that in the continuum limit the gauge condition (p3D corresponds 
to a back ground gauge 

[d, + iAl{x),A^{x)] = 0, (35) 

where the back ground gauge field A^[x) is an optimally chosen center vortex 
field (see [|^ for details). In the absence of center vortices in the considered 
gauge field Afj_{x), the gauge (^) coincides with the Landau gauge. 

The gluon form factor in MCG is shown in figure ^ By comparing figures |^ 
and ^ we observe that the form factors in Landau gauge and MCG drastically 
differ in the intermediate region. This is because in Landau gauge most of the 
information is accumulated in the adjoint (coset) part of the links while in the 
MCG the part of the non-perturbative content is shifted to its center part. 
Furthermore, comparing figures and ^ we find that the gluon propagator 
in MCG basically agrees with the one where the center vortices are removed 
and the Landau gauge is subsequently implemented. This underlines the 
observation that in MCG confinement physics is relegated to center vortices. 

3.5 The Gribov noise 

Finally, let us check how strongly the gluonic form factor -F(p^) depends on 
the choice of gauge, i.e. on the sample of maxima of the variational condition 
([Tl|) selected by the algorithm. For this purpose, we adopt an extreme point 
of view by comparing the gauge implemented by the iteration over-relaxation 
(10) algorithm with the gauge obtained by simulated annealing (SA). The 
results obtained for the form factor in both cases are shown in figure ^. We 
find, in agreement with pSf, that, in the case of the gluonic form factor, the 
Gribov noise is comparable with the statistical noise for data generated with 
(5 e [2.1,2.5] (scaling window). 

4 Gluonic spectral functions 
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Figure 4: Tlie gluonic form factor -F(p^) for a 10^ x 20 lattice in the gauge 
10 and SA, respectively (left panel) and compared with previous results (10, 
16^ X 32) (right panel). 

4.1 The spectral density 

Let us consider the following spectral representation of the trace of the Eu- 
clidean propagator DijP') 

pirn) 



D{p' 



dm ■ 



+ m^ 



(36) 



'0 P 

where p{m) is the spectral density. For example, the propagator of free par- 
ticle with mass nip is represented by a spectral density p{m) = 5{m — nip). 
In general, the spectral density carries information on the strength with 
which single particle states contribute to the correlation function of interest. 
Hence, the spectral density must be positive if the space of physical parti- 
cles is considered. Note, however, that this constraint must be abandoned 
if un-physical (while gauge dependent) correlation functions, e.g. the gluon 
propagator, is investigated. Moreover, the so-called negative norm states 
play an important role in Yang-Mills theory to circumvent the cluster de- 
composition theorem, what is inevitable to accomplish confinement of the 
theory. For a more detailed discussion of these issues see |T^ . 



In order to study the contribution of hypothetical negative norm states to 
the gluon propagator, we will calculate the spectral density from the gluonic 
form factor 



F{p' 



dmp{m) 



p^ 

p2 _|_ jy^2 



(37) 
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Taking the derivative of F(p^) ( [T6| ) with respect to p, one finds 



m 



2 



F\v^) = 2p / dmpim) "" . (3^ 

J [p^ + m^) 

Hence, one immediately concludes that the form factor would be a monotonic 
function of the momentum if only positive norm states contribute, i.e. if 
p(m) > 0. An inspection of the results of the previous section shows that 
this is not the case for the gluon propagator. 

Further insights are provided by the sum rules which are obtained from a 
large momentum expansion of the form factor (0), i.e. 

oo 

F(/) = Y, ^n^ ' ^n = (-1)" / dm p{m) m^" . (39) 

where n < rtc'is restricted from above in order to guarantee the existence of 
the momentum integrals. In particular, we observe that 

lim -F(p^) = / dm p{m) . (40) 

p2^00 J 

There is is a substantial difference between the spectral function of a free 
massive particle and the one of Yang- Mills theory: while in the free particle 
case the form factor approaches unity at large momentum transfer, we know 
from perturbation theory that the Yang-Mills form factor vanishes at large 
momentum transfer. In view of (^0[) , the crucial observation is that the free 
particle form factor F(p^ — i> oo) = 1 is consistent with a positive definite 
spectral function, while the vanishing form factor F{p^ — »• cxo) = in the 
Yang-Mills case implies that unless p{m) vanishes identically it must change 
sign. 

The remaining two subsections are devoted to numerically estimate the spec- 
tral function p{m) from the form factor data in oder to get insights into the 
qualitative behavior of the spectral function of a confining theory. 

4.2 The generalized Maximal Entropy Method 

In order to extract the spectral function p{m) from our gluon propagator 
obtained in the lattice simulations we employ the Maximal Entropy Method 
(MEM). There is a wide span of applications for the MEM ranging from image 
modeling to solid state physics [Q. It was recently applied to reconstruct 



the spectral density of mesons from correlation functions obtained in lattice 
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calculations |^. Since these applications use the positivity of the spectral 
density, a slight generalization of the MEM is necessary for our purposes. 



Given the definition 



fuEMiP 



dm p{m] 



p" 



p2 _|_ ^2 



we define the MEM potential functional by 



V\p] 



dp 



a{p) l 



-^(P^) - fuEMiP^ 



(41) 



(42) 



where a{p) denotes the standard deviations of the measured values FijP'). 
For given F{p'^)^ minimization of the functional V[p\ with respect to the 
spectral function p{m) defined by ( ^T]) corresponds to a least square fit and 
would result in the optimal choice for p{m) representing the data. However, it 
appears that substantial changes of the function p{m) produce minor changes 
of the function fuEMip^) (PD which are comparable in size with the error 
bars of the measured function FijP'). In practice, the numerical algorithm 
which tries to minimize the potential (021) is unstable. 



(43) 



For circumventing this problem, one defines a MEM action functional 



S[P\ 



" '-'entropy [PJ 



+ V[p\ 



where the parameter a regulates the influence of the entropy function Sentropy[p\ 
(to be specified below) on the minimum of the functional S[p]. The action 
functional (^) is minimized instead of the potential (^21) . Thereby, the en- 
tropy factor S^ntropyip] IS a measure for the deviation of the spectral density 
p{m) from a default density pde/(^7i). Whenever the potential term (^) is 
not conclusive on the precise form of the spectral density due to the size 
of the error bars, the numerical algorithm minimizes the deviations of p{m) 
from the default model specified by pdefijn) ■ The entropy factor S^ntropylp] 
is also used to encode constraints on the function p{m). For instance, if one 
would like to insist on the positivity of the spectral function, the generic 
choice is [P] 



S, 



entropy 



[P] 



dm 



p{m) log 



p{m) 



since the functional 



pdefKm) 
has a unique minimum for p{m) 



p{m) 



(44) 



Pdef[fn) 



For each value a a minimum of the action (|43|) exists for properly chosen en- 
tropy functionals. Let paijn) denote the function which minimizes (|45|), and 



20 



let s{a) be the corresponding minimuni value. The most probable spectral 
function is then defined by 



Pav[m) 



—. I da PaiTTi) exp< —s{a) > , Af := da exp< —s{a) > . 

(45) 
Of course, it is mandatory to check that fuEMip^) (HI) which is calculated 
with Pav{fn) does indeed represent the data F{p^) within statistical error 
bars. 
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Figure 5: The spectral functions obtained with the entropy functionals S**^^-* 
and S^'^\ respectively, (left panel) and the corresponding form factors (right 
panel) . 



Since we must abandon the positivity constraint, we cannot use the standard 
entropy functional (^4]). We will here explore the two functionals 



9^^^ \o\ 

'-^entropy Lr J 



s. 



(4) 
entropy 



[P] 



dm 
dm 



dp{m) 
dm, 



-\ 2 



p{m) - pdef{m) 



n2 



(46) 
(47) 



where we use a "default model" which suppresses the density p{m,) at small 
and at high values of m, e.g. 



Pdef[rn) 



m 



m^ + (2 rrip) 



(4^ 



where the mass scale is chosen to be m,r 



1 GeV. The first functional 



suppresses large gradients and, hence, ensures that the spectral function is 
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Figure 6: Same as figure |^ for tlie case of the non-confining Yang-Mills 
ensemble. 



smooth. The second functional allows for large gradients, but minimizes 
the deviation of the spectral function from that of a "default model". The 
two functionals have been tested for the case of the form factor of a free 
particle (see appendix 0). The spectral function of the cases ( ^Tj) and (^51), 
respectively, have been compared with the result produced with the standard 
entropy functional (|4^) . Since less information is supplied for the choices (0, 
481) , the exact infinite volume density is reproduced to less accuracy than in 
the case using the standard entropy functional (PI). These results stress the 
importance of the physical constraints on the spectral function. We point 
out that, unless additional information on the spectral function of the gluon 
form factor is obtained and is used to constrain the latter, one must accept 
the uncertainty which becomes visible if two different entropy functionals 
(such as ( ^TD and (|i8| ) are used. 



4.3 MEM fit of the lattice form factor 

In a first step, we do not attempt to supplement additional physical informa- 
tion to the MEM approach, but produce a direct MEM fit of the form factor 
data by using the two completely different entropy functionals (^Uj), (|i7| ) in 
order to get a clue about the residual freedom in the choice of the spectral 
function p{m). Our findings for the spectral function corresponding to the 
form factor F{q^) of the gluon propagator in Landau gauge is summarized in 
figure ^ The MEM approaches which employ the entropy functionals (^61), 
(|T7|) produce qualitatively the same result for the spectral function: p{m) is 
a rapidly oscillating function of m. 
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Figure 7: Same as figure ^ for tlie case tliat tlie perturbative tail lias been 
added (full Yang-Mills ensemble) . 

We compare these findings with the ones obtained from the modified Yang- 
Mills ensemble from which the center vortices have been removed by the 
method described in the previous section. As shown there, the modified en- 
semble is non-confining. Since in this case, as in the case of the full ensemble, 
the form factor is non-monotonic due to the perturbative tail, negative norm 
states are inevitable. The MEM fit to the data indeed shows an oscillation 
of the spectral function p{m) the amplitude of which is, however, orders of 
magnitude smaller than in the case of the full Yang-Mills theory (see figure ^. 

Finally, we will assume that the perturbative regime is already approached 
at the upper limit of the momentum range which has been explored by the 
lattice calculation so far. Therefore, we will artificially extend the momen- 
tum range of the form factor by attaching the perturbative tail to the high 
momentum form factor obtained from the lattice data. In addition, we add 
a Gaussian noise to the perturbative tail. The noise is of the same order of 
magnitude as the statistical error bars, and specifies the degree of freedom 
which the MEM method might exploit to reproduce the form factor data. 

Figure |^ shows the result for the spectral function for the case of the full 
Yang-Mills gluonic form factor. Most important is the observations that 
the oscillations are stable in the position. There is a enhancement of the 
amplitudes which is naturally expected if the momentum range of the data 
grows. 

We also repeated the analysis of the spectral function for the case of the 
non-confining ensemble (see figure §). Also in this case, we find that the 



23 



B 







A 




— S'^dat 


30 


\i\ 






- S'^'.dat 

S'", dat+per 

S'"', dat+per 


10 
10 


1 

1 


\ "^ / 


^^-~- 


30 


V' 

y 






- 













perturbative 
tail 


a. 


1.5 

1 




^ 


% 


^--^ 




0.5 








o dat+per 

g(.. 

s'"' 



2 



Figure 8: Same as figure ^ for tlie case tliat tlie perturbative tail lias been 
supplemented to the lattice data (modified non-confining ensemble). 

attachment of the known perturbative behavior at large momentum transfer 
does not qualitatively change the behavior of the spectral function p{m), but 
results in an enhancement of the amplitude. 

4.4 Infra-Red enhancement and negative norm states 
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Figure 9: Spectral functions corresponding to the model form factors Fi^ijP') 
1) for several values k of the infra-red strength. 



In order to detect the origin of the rapid oscillations in the spectral density 
pirn) for the case of the full gluonic from factor (see figure ||), we are guided 



24 



by the mass fit (^) of tlie lattice study and take tlie class of functions, i.e 

.2 . ..,,4 1 



F.{p' 



p- 



p^ + mi 



KjJ 



p^ + m\ 



+ 



^09V^^'y. 



13/22 



(49) 



where mi = 0.64 /z, m^ = rni = /^ (varying k), as a model for the gluonic form 
factor. Here, k controls the medium infra-red strength of the form factor. 
The MEM of the form factor fits which employ the entropy functional S''-^^ 
are shown in figure ^ for k = 1, 0.33, 0.11. These fits of the form factor fall on 
top of model functions -Fk(p^) within plotting accuracy. The corresponding 
spectral functions are also shown in p. One observes that amplitude of the 
oscillations of p{m) increases if the infra-red strength, i.e. k, rises. 

5 Conclusions 

We have studied the gluonic form factor F{p^) in Landau gauge and maximal 
center gauge, respectively, by means of the continuum extrapolated SU(2) 
lattice gauge theory. The gauge field has been defined from the adjoint link 
variables. It is well known [^] that the Landau gauge condition (|llD does 
not completely fix the gauge, but one must select a subset of configurations 
lying within the first Gribov horizon. We have studied two such possibilities: 
firstly, an iteration over-relaxation algorithm (10 gauge) which averages over 
the configurations within the first Gribov horizon, and, secondly, simulated 
annealing which selects the configuration belonging to the fundamental mod- 
ular region (SA gauge). Within the statistical error bars, both methods have 
yielded the same gluonic form factor. Our results are also in agreement with 



the findings using the Laplacian Landau gauge [^ where the problem of 
Gribov ambiguities is evaded in an elegant way. 

Close to zero momentum transfer, the gluonic form factor is mass dominated 
where the mass is ~ 650ib20 MeV (we used \/a = 440 MeV as reference). The 
uncertainty is due to statistical errors as well as due to Gribov ambiguities, 
i.e. using 10 and SA-gauge, respectively. Our numerical results are stable 
against a variation of the UV-cutoff (i.e.vr/a) and the physical volume by 
roughly an order of magnitude (see table |I|). The data points of the form 
factor in the far infra-red were obtained by using lattice sizes of L ?a 9 fm. 
This volume corresponds to a topological mass 27r/L ^ 140 MeV which is 
much smaller than the dynamical gluon mass mi ^ 650 MeV. This explains 
the small finite size effects on the form factor observed by our numerical 
simulation. 
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In the medium momentum range we have found a rather pronounced peak 
while at high momenta our numerical data nicely reproduce the result from 
perturbative Yang-Mills theory. Our numerical data are well fitted over the 
whole momentum range by the formula (p5|) which might be useful for further 
phenomenological oriented investigations. 

Our SU(2) gluonic form factor in Landau gauge shows the same qualita- 
tive behavior than its SU(3) counterpart using Landau gauge |^ or using 
Laplacian Landau gauge ||35|, ^ . 

A focal point of our studies is the information on quark confinement which 
might be encoded in the gluon propagator in Landau gauge. By removing the 
confining vortices from the ensemble by hand, we are left with an ensemble 
which does not confine quarks (see figure ^ left panel). After the implemen- 
tation of Landau gauge, we have seen that a good deal of strength is removed 
in the medium momentum range. We have therefore established a relation of 
the infra-red strength of the gluonic form factor in Landau gauge and quark 
confinement. 

We have compared our result for the gluonic form factor with that obtained 
by solving a truncated set of Dyson-Schwinger equations (DSE) in continuum 



formulated Yang-Mills theory |31|. The DSE result is in qualitative agree- 
ment with our findings, but does not reproduce the data on a quantitative 
level (except in high momentum, perturbative regime): the power law be- 
havior of the form factor in the vicinity of zero momentum is different, and 
the peak of the intermediate momentum range is less pronounced given that 
both approaches give same results in perturbative momentum regime. 

Finally, we have studied the gluonic spectral functions which reproduce the 
numerical data for the gluon form factor. Using sum rule techniques, it 
becomes clear that the spectral function necessarily comprises negative parts. 
Using a generalized version of the Maximum Entropy Method (MEM), we 
have found that this is the case for the full form factor as well as for the form 
factor obtained for the non-confining ensemble. At a quantitative level, the 
spectral function corresponding to the full form factor shows large amplitude 
fluctuations while the amplitude of the oscillations in the case of the spectral 
function obtained from the non-confining ensemble are moderate. 
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A Gauge fixing on the lattice 

In practical calculations, the following convenient procedure, which circum- 
vents the explicit evaluation of the Faddeev Popov determinant, is adopted: 
a number of {Ufj,{x)}i=i,_,n of statistically independent ensembles are gener- 
ated by standard techniques. Each of these ensembles is then subject of the 
implementation of the gauge condition (p^, i.e. 



{U,{x)h - {U^{x)h , Vz . (50) 

Depending on the algorithm, a single candidate of all possibile maxima of 
the variational gauge condition is randomly selected (see discussion in subsec- 
tion |2l^ ). Independently of the procedure which selects the gauge transfor- 
mation Q{x), a particular member {U^{x)}k, /c G {1 . . . n} of the gauge fixed 
configurations is generated with a frequency proportional to the probability 
distribution of the gauge fixed sub-manifold. An estimator of a quantity A 
is obtained by 

A^\^iVUl . (51) 

Let us consider the particular example that A is a gauge invariant combina- 
tion of fields, A = A^. One obtains with (Ell) 



n n n 

^)^-EM = -EM = -EM -(52) 

I I I ^ 

Hence, gauge invariant quantities which are calculated from the gauged con- 
figurations evidently coincide with the ones obtained from un-fixed configu- 
rations. 



B Generalized MEM method (example) 

We will briefly investigate several proposals for the MEM entropy functional 
Sentropy[p\ VO- (0) for the case of the form factor of a free particle of mass nip 
where noise has been supplemented to the form factor by hand to simulate 
statistical uncertainties. The form factor is given by 

np') = ^^ ■ (53) 

In order to illustrate the sensitivity of the spectral function p{m) to the con- 
straints, which are incorporated in Sentropy[p], we will explore four different 
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least square fit (case S 




Figure 10: Possible spectral functions (differing by constraints) which fit 
the form factor of a free massive particle. 



entropy functionals, i.e. 

9^^) \o\ - 

^entropy iF\ 



s. 



s, 



(2) 
entropy 

(3) 
entropy 



[P] 
[P] 



9^"^) In] 

'^entropy Irl 



dm 
dm 
dm, 



dpim) 
. dm 

p[m) log — - — p[m) 



p{m) 



Pdef['m) 
Pdefifn) 



(54) 
(55) 

(56) 

(57) 



In order to perform the MEM program outlined in subsection [4.2| , all integrals 
are converted Riemann sums. We used 135 data points representing F(p^) 
and Nm points to represent the spectral function p{m,). 



The choice S, 



(1) 

entropy 



corresponds to a least square fit of the spectral sum 



fuEMip^) to the "measured" data. In this case, it turns out that the set 
of points p{m) which represent the data set best is not a smooth function 
(see figure |10| for the case m^ = 1 and A^^ = 90). 

The minimal requirement which we wish to incorporate in the spectral func- 
tion p{m) is that p{m) is a smooth function of m. This is achieved be the 



entropy functional S\ 



(1) 

entropy 



p5| ) which disfavors functions with large gradi- 



ents. The spectral function is also shown in figure |10] where we have chosen 
Nm = 90. One observes that a peak develops at m = rup approximating the 
exact infinite volume and zero noise spectral function pexact = ^i^n — m^). 
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Figure 11: The form factor of a free massive particle fitted by several spectral 
functions. 

Let us compare, this result with the result of the standard approach de- 
manding positivity. For this purpose, we use SlJ^^^y (0). This choice for 
the entropy functional does not constrain gradients. The smoothness oi p{m) 
is here incorporated by minimizing the difference to a smooth default density 
Pdefijn)- Here, we used 



Pdef{m) 



m 



m^+ {2mpy 



(58) 



which ensures that the spectral function vanishes for m — > and m -^ 
cxD, respectively. This is the standard MEM approach to form factors with 
positive definite spectral functions. 

Finally, we use an "Euclidean" norm to measure the difference of the spectral 



function with the default model (see S, 



(4) 
entropy 



[P] 



Thus, the constraint 



to positivity is abandon, and we, hence, expect that agreement of the MEM 
estimate of the spectral function with the exact result becomes worse. This 
is indeed observed (see figure ITUI). 



We point out that all above MEM suggestions for the spectral function p{m) 
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fit the form factor F{p'^) witliin the (artificial) statistical error bars. The 
best result for the spectral function is achieved with the entropy functional 
^entropyip] (0) which incorporates the most severe constraints, i.e. smooth- 
ness and positivity. 
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